The polynomial invariants of four qubits 



Jean-Gabriel Luqu^ and Jean- Yves ThiborOl 
Institut Gaspard Monge, Universite de Marne-la-Vallee 
77454 Marne-la-Vallee cedex, France 
(Dated: December 14, 2002) 

We describe explicitly the algebra of polynomial functions on the Hilbert space of four qubit states 
which are invariant under the SLOCC group 5'L(2,C)*. From this description, we obtain a closed 
formula for the hyperdeterminant in terms of low degree invariants. 

PACS numbers: 03.67.Hk, 03.65.Ud, 03.65.Fd 



I. INTRODUCTION 

Various classifications of states with up to four qubits 
have been recently proposed, with the aim of understand- 
ing the different ways in which multipartite systems can 
be entangled 0, 0, 0, 0, However, one cannot ex- 
pect that such classifications will be worked out for an 
arbitrary number k of qubits, and there is a need for a 
coarser classification scheme which would be computable 
for general k. In 6], Klyachko proposed to assimilate en- 
tanglement with the notion of semi-stability of geometric 
invariant theory. In this context, a semi-stable state is 
one which can be separated from by a polynomial in- 
variant of SL{2, C )^ , the point in the geometric approach 
being that explicit knowledge of the invariants is in prin- 
ciple not necessary to check this property. 

In this paper, we construct a complete set of algebraic 
invariants of 4-qubit states. This allows us to identify 
the semi-stable states in the classification of Verstraete 
et al. 3, and to obtain a simple closed form for the 
hyperdeterminant . 

Let = be the local Hilbert space of a spin i 
particle, and Ti — be the state space of four parti- 
cles, regarded as the natural representation of the group 
G = SL{2, C)^, known in the context of quantum infor- 
mation theory (QIT) as the group of reversible stochastic 
local quantum operations assisted by classical communi- 
cation (SLOCC) [10. 

If b); j = 0,1 is any basis of V, a state l^*) can be 
written as 



^vki\i) li) \k) (S> \l) 



(1) 



and the question of which normal form can be achieved 
for I^P) by varying independently the bases of the four 
copies of V has been solved only recently 3 , although 
the case of three-qubit states is classical and relatively 
simple 

In the following, we give a complete description 
of the polynomial functions f{Aijki) which are invariant 
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under the SLOCC group SL(2, C This amounts to the 
construction of a moduli space for four qubit states. Our 
strategy is to find first the Hilbert series of the algebra 
of invariants ^7. Next, we construct by classical methods 
four invariants of the required degrees. The knowledge 
of the Hilbert series reduces then the proof of algebraic 
independence and completeness to simple verifications. 
The values of the invariants on the orbits of 0| are tab- 
ulated in the Appendix. 



II. THE HILBERT SERIES 

Let J7d be the space of G-invariant homogeneous poly- 
nomial functions of degree d in the variables Aijki ■ Using 
some elementary representation theory, it is not difficult 
to show that JT^ is zero for d odd, and that for d = 2m 
even, the dimension of is equal to the multiplicity of 
the trivial character of the symmetric group &2m in the 
fourth power (x™™)'' of its irreducible character corre- 
sponding to the partition [m, m]. This is the same as the 
scalar product ((x'"'")^|(x'"'")i)jwhich can be evaluated 
by means of the formulas of H, Q giving the decompo- 
sition into irreducible characters of any product x^X^ 
when A and ^ have at most two parts. This yields the 
Hilbert series of J" = 0^ in the form 



d>0 



(l-t2)(l_i4)2(i_^6) • 



(2) 



This formula shows that Conjecture 2.6.5.3 of |Q| cannot 
be correct, since it predicts that the hyperdeterminant, 
which is of degree 24, should be one of the generators. 
Actually, the algebra of invariants is free on generators 
of degrees 2, 4, 4, 6, as suggested by the Hilbert series. 



III. A FUNDAMENTAL SET OF INVARIANTS 

Indeed, it is possible to construct invariants of the re- 
quired degrees and to check that they are algebraically 
independent. To reduce the size of the expressions, we 
shall write the components of \^) as 



A 



ijkl — ■ 



0, 



15, 



(3) 



2 



where r is the integer whose binary expression is ijfc/, 
that is, r = 8i + Aj + 2k + I. We shah consider them as 
the coefficients of a quadriUnear form 

1 

A(x,y, z,t)= ^ AijuXiVjZkU (4) 

onVxVxVxV. Such a form is known to have an 
invariant H of degree 2, which is also one of the hyper- 
determinants introduced by Cayley [2J|. It is given by 

iJ = floais - aiai4 - 02013 + 03012 

-04011 + 05010 + OgOg - O7O8 , (5) 

and the two independent invariants of degree 4 are any 
two of the 3 determinants which can be formed by inter- 
preting A as a hnear map ^ C (see [13) 
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One has the relation 



Oq 


04 


Og 


ai2 


Ol 


05 


09 


ai3 


0.2 


06 


aio 


Oi4 


03 


07 


On 


ai5 


Oq 


Os 


a2 


aw 


Ol 


09 


as 


On 


O4 


ai2 


Og 


Oi4 


05 


ai3 


07 


ai5 


Oo 


Ol 


Og 


09 


02 


03 




On 


04 


05 


ai2 


ai3 


06 


07 


Oi4 


ai5 



L + M + N = 0, 



(6) 

(7) 

(8) 
(9) 

(10) 



but it is easily checked that any two of them are linearly 
independent, and also that cannot be expressed as a 
linear combination of them. 

To construct a sextic invariant algebraically indepen- 
dent from the previous ones, we shall apply the methods 
of classical invariant theory, and first find some covari- 
ants, that is, homogeneous G-invariant polynomials in 
the form coefficients Aijki and in the original variables 
(see, e.g., for a modern presentation). The dimen- 
sion of the space Cd,fci,/c2,fe3,fe4 of covariants which are of 
degree d in A, ki in x, and so on, is equal to the multi- 
plicity of the trivial character of 6 ^ in the product 



^h+ki ,li ^l2+k2d2 ^l3+k3,la ^h+ki,li 



(11) 



where d = 2li + ki for all i. This can still be evaluated 
from the knowledge of the products x^x'^ of two-part 
characters, and one can see that A has six covariants of 
degree 2, which are biquadratic forms in all possible pairs 
of variables. Such covariants are easily constructed, these 
are the determinants of order 2 of the partial derivatives 



of A with respect to the complementary pair of variables, 
e.g., 



6a;y(x,y) = det 



d^A 



(12) 



Each of these biquadratic forms can be interpreted as 
a bilinear form on the three-dimensional space 5*^(0^), 
and one can define 3x3 matrices by 



yovi 
vl 



(13) 



and similarly for the other pairs. Let, for any pair of 
vector variables (u, v). 



= det(i3™) . 



(14) 



These determinants are sextic invariants of A. Since the 
space of sextic invariants is four-dimensional, they must 
be linearly dependent. In fact. 



Dyt and D^t 



(15) 



but Dxy, Dxz and Dxt are linearly independent. One can 
check that 



HL ^D,, 
HM ^D,t' 

HN =D,y 



(16) 
(17) 
(18) 



and that is not in the subspace spanned by the D^v 
The above results, together with the knowledge of the 
Hilbert series, allows now to prove that the algebra of 
invariants is free, and that any of the D^v^s can be taken 
as the generator of degree 6. Indeed, it is sufficient to 
check that the Jacobian matrix of the choosen generators 
has rank 4 (this is easily done using the specialization 
Gated of the Appendix). 
We will use in the sequel 



J^C[H, L, M, D^t] 



(19) 



IV. THE HYPERDETERMINANT IN TERMS 
OF THE FUNDAMENTAL INVARIANTS 

Here, according to the general formula of [l2| , the Cay- 
ley hyperdeterminant (in the sense of |E^3) is of degree 
24. It must therefore admit an expression in terms of the 
fundamental invariants, whose explicitation is an inter- 
esting question. To answer it, we shall need again the 
covariants buv Let us use, for example, bxt- We can con- 
sider ^ as a trilinear form T in x, y, z, the fourth variable 
t being treated as a parameter. The Cayley hyperdeter- 
minant Det(T) of this trilinear form is homogeneous of 
degree 4 in its coefficients, which are themselves linear 
forms in t. Hence, i?(t) = Det(T) is a binary quartic 
in to,ti, and we can form its discriminant A, which will 



3 



be an invariant of A. According to Schlafli [13 (see also 
[E^|)i in this case, A is equal to Det{A). 

It follows from the well-understood invariant theory of 
binary trilinear forms 13, Ts'l that i?(t) is equal to the 
discriminant of the quadratic form in x 



Qt(x) = 6j:t(x,t) , 



that is 



Let 



i?(t) = det 



\dxidxj 



(20) 



(21) 



R{t) = coto + 4ciigti -I- &C2tltl + Uotl + Cit\ . (22) 

It is well-known that the algebra of invariants of the bi- 
nary quartic is free over the two generators 



2 ; 



S = C0C4 — 4ciC3 + 3c: 

T — C0C2C4 - C0C3 -I- 2C1C2C3 - C^C4 - 
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and that its discriminant is given by 
A^S^-27T^. 



(23) 
(24) 

(25) 



In the classical language, S is the apolar of R with itself, 
and T is its catalecticant (see [lljl. 

The invariants S and T of i? being obviously invari- 
ants of A, the problem of expressing Det (A) is terms of 
the fundamental invariants of A is reduced to the one of 
finding the expressions of S and T [25j| . 

With the help of a computer algebra system, we obtain 
the values 



S = —H^--H^L+-H^M-2HD^t 
12 3 3 

-^^(L^ + LA/ + M2) 
o 



(26) 



and 



216 18 ^ J Q xt 



LM + 2M 



/r2^ 



-H{L - M)D,t 



-LM{L-M)+Dl, 



(27) 



Setting D = D^t, U = ~ 4:{L ~ M) and V = 
12{HD — 2LM), these expressions can be recast into the 
more elegant form 



125 = ~2V , 
216r = - mV + 2161)2 . 



(28) 
(29) 



This suggests that U and V might have a geometric 
meaning. Actually, similar expressions occur in the 
course of Schlafli's calculations [lj|. He does not mention 
their invariant theoretic meaning, however, and prefers 
to end up with an expression of A as a polynomial in 
H, W = D,y + D,, + D,t, ^ = + AP + and 
11= [L — M) {M — N) {N — L) , which are invariant under 
permutations of the indices ijkl. 



V. CONCLUSION 

A fundamental issue in QIT is the understanding of 
entanglement. However, as pointed out in Q, there is no 
universal agreement on the precise definition of entan- 
glement and on what should be its proper measure. It is 
apparently this lacune which motivated recent attempts 
to obtain a complete classification of fc-qubit states under 
the SLOCC group G [JIIH. 

Some familiarity with classical invariant theory leaves 
little hope that such a classification can be achieved in 
general. If we compare with the somewhat easier clas- 
sical problem of binary forms, which, in physical terms, 
amounts to the classification of single spin s states un- 
der 5i(2,C), a complete solution is known only up to 
spin s = 4 (with still some unsolved questions in the case 
s = 7/2), and most experts agree that the other cases 
will remain out of reach. 

So, it is unlikely that the complete SLOCC classifica- 
tion of fc-qubit states will ever be found for fc > 8, and 
it is probable that formidable c omp utational difficulties 
will arise well before this value j2y|. Actually, the orbit 
structure is still completely unknown for fc > 4. 

Now, if we adopt the definition of entanglement pro- 
posed in 0, that is, to identify entangled states with 
the semi-stable vectors of geometric invariant theory, the 
main result of the present paper can be interpreted as 
a numerical criterion of entanglement for 4-qubit states. 
Indeed, a semi-stable state is by definition a state which 
can be separated from by some invariant polynomial. 
Thus, according to 0, an entangled 4-qubit state would 
be one for which at least one of the four polynomials 
H, L, M, D takes a nonzero value. As we shall see be- 
low, this definition needs to be improved. However, it 
is plausible that refined entanglement measures for four- 
qubit states might be built from the absolute values of 
these invariants. These would be natural generalizations 
of the concurrence G and the 3-tangle r in the two and 
three qubit cases, which are proportional repectively to 
the absolute values of the determinant and of the hyper- 
determinant 0, the only polynomial invariants in these 
cases. 

From a geometric point of view, our results show that 
the moduli space of entangled states is the weighted pro- 
jective space P(l,2,2,3), which can be embedded as a 
rational threefold in 13-dimensional projective space. Of 
course, the approach to semi-stability and moduli spaces 
by explicit construction of the polynomial invariants has 
its limits, and it is unlikely that this can be done for more 
that 5 qubits. 

Note also that the notion of semi-stability can be used 
only to characterize some generic kind of entanglement. 
Indeed, even in the three-qubit cases, the so-called W- 
state ^(|001) + |010) + |100)) is not semi-stable (its hy- 
perdeterminant is zero), although it should certainly be 
considered as entangled (even in a strong sense, accord- 
ing to U). The natural candidates for constructing fur- 
ther measures of entanglement appropriate to such states 
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are the covariants of classical invariant the ory , which are 
completely known in the three-qubit case [I3 • 

We expect to be able to describe in a forthcoming pa- 
per the algebra of covariants in the 4-qubit case, which 
would not only reproduce the complete classification of 
the orbits, but also to give the equations of their clo- 
sures, which are algebraic varieties, and provide new in- 
sights about entanglement measures for unstable states. 
Another case whichs seems to be readily tractable is the 
case of three spin 1 states. The geometric classification 
of the orbits was known by 1938 [Ig, and numerical cal- 
culation of the Hilbcrt series of invariants up to degree 
108 indicates that it should be 



and the hyperdeterminant is -^V {a'^ , b'^ , , (P)"^ , where 



h{t) 



1 



(1 -t6)(l _f9)(l _tl2) 



(30) 



Since independent invariants Ig, Ig, I12 of the appropriate 
degrees are known 0, one can consider that the SLOCC 
classification and the semi-stability problems are essen- 
tially solved in this case. 

Other cases of less practical importance, such as those 
including two spin ^ particles and one particle of spin 
s > 1, are easily solved. For s = 1, there is only one 
invariant of degree 6, the hyperdeterminant. For s — 3/2, 
the hyperdeterminant is identically zero, but there is still 
one invariant of degree 4, which is the only determinant 
that can be formed by displaying the components of \'^) 
in a 4 X 4 matrix. Finally, for s > 2, there are no invariant 
polynomials. 
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APPENDIX A: APPLICATION TO THE 
CLASSIFICATION OF VERSTRAETE ET AL. 

To conclude, let us discuss the semi-stability of the 
orbits obtained in Q (see this reference for notation) . For 
the family Gabcd, the values of the fundamental invariants 
are HI 



and A = 0. For L„ 
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and A = 0. For -^0203^1 j there is still one nonzero invari- 
ant 
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D = ~-{ad-bc){ac~bd){ab-cd) (A4) 



H = a\ 



(A21) 



so that all the above 6 families of orbits are semi-stable, 
whilst the remaining three are unstable. 
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